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sampled linearly, which are found to be structurally critical. However, we encounter 
interesting differences both in the theory (where the role of the history frequency 
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of the models. Our theoretical predictions are supported by numerical simulations. 
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1. Introduction 

Minority Games (MG) [U [2] are relatively simple and transparent models that were designed 
to increase our understanding of the complex collective processes which result from inductive 
decision making by large numbers of interacting agents in simplified 'markets'. They can be 
seen as mathematical implementations of the so-called El Parol bar problem [3] . Many versions 
of the MG have by now been studied in the literature. They differ in the type of microscopic 
dynamics used (e.g. batch versus on-line, stochastic versus deterministic), in the definition of 
the information provided to the agents (real-valued versus discrete, true versus fake market 
history) and the agents' decision making strategies, and also in the specific recipe defined for 
the conversion of the external information into a trading action (inner products versus look-up 
tables). In MG models with 'fake' market memory, proposed first in [4], at each point in time 
all agents are given random data upon which to base their decisions, rather than the actual 
market history. They have the advantage of being Markovian and were therefore the first to be 
studied and solved in the theoretical physics literature using techniques from equilibrium E] 
and non-equilibrium [3 [U |9] statistical mechanics, as first developed for and applied to spin- 
glasses. Models with true market history are non-Markovian and therefore more demanding, 
but also here theoretical progress is being made |10^ lllj. For a more extensive introduction to 
MGs and their statistical mechanical theory we refer to the recent textbooks \12\ [T3] . 

To our knowledge, however, there has not yet been any such exact solution for MG 
models in which the decision making strategies are defined via inner products rather than 
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look-up tables |14t I15j: not even for the simplest Markovian case, where the market histories 
are fake. Inner product MGs have so far tended to be studied either numerically, or at the 
level of investigating properties of the stochastic microscopic laws, in spite of the fact that 
one might well argue that the definitions of the 'inner product MGs' are possibly closer to 
how financial time series tend to be predicted by practitioners |19j (as these usually involve 
generalized linear market prediction models). Moreover, in fact not a single study has so far 
been published of inner product MGs with real as opposed to fake histories, there are not even 
papers based on numerical simulations only. 

In this paper we first generalize the existing generating functional theory [TI] that was 
developed for look-up table MGs with real and/or fake histories to a larger class of MG 
versions, which includes the familiar look-up table MGs and inner product MGs as special 
cases. We then focus on the analysis of the dynamics and phase transitions of inner product 
MG versions with and without real histories, and show that in the infinite system size limit: 
(i) for inconsistent fake histories the inner product and lookup table MG versions behave fully 
identically, including short time correlation and response functions and the volatility (apart 
from a trivial transformation related to time-scales), (ii) similar to what was found in [11] 
for look-up table MGs, for inner product MGs the introduction of real histories affects the 
observables in the ergodic region, but does not generally change the phase diagram, except 
when the histories are sampled by unbounded functions (see below), (iii) the key observable 
to measure the differences between real and fake histories in our generalized MG models is 
the spectrum of the history covariance matrix in the stationary state, which for inner product 
MGs reduces to the history frequency distribution, (iv) for inner product MGs an approximate 
calculation of this spectrum leads to a fully closed stationary state theory, the predictions of 
which are supported by numerical simulations, and (v) inner product MGs with predominantly 
real histories where the inner product is taken between strategy entries and the actual values 
of the past market bids (or, more generally, by an unbounded function of these bids rather 
than e.g. a saturating function) are permanently critical, giving a volatility that is either zero 
or infinity (separated by a novel critical value for the control parameter a, unrelated to the 
standard ergodicity breaking transition point). 

2. Definitions 

We imagine having a system of A'" agents, labeled hy i = 1, . . . , N. At each step £ S {0, 1,2,.. .} 
of the game each agent submits a 'bid' bi{i) S IR to the market. The (re-scaled) cumulative 
bid at stage i is defined as 

A{£) = ^Y.b^ie) + A,{e) (1) 

Here Ae{i) denotes an external contribution that represents e.g. random market perturbations 
or actions by market regulators. Each agent i determines his bid bi{£) at each step £ on the 
basis of external information, which is to represent historic market data, using his preferred 
decision making strategy at that point in time. Each agent has S such strategies, labeled by 
o = 1, . . . , S*. Profit is made by those agents who find themselves subsequently in the minority 
group, i.e. in the case of A(£) > by those with bi{£) < and in the case of A{i) < by those 
with bi{£) > 0. 

MG versions differ in their definition of the information provided to the agents and in 
how this information is converted into trading decisions. In the original MG, the information 
consisted of the signs of the overall market bids A(£) over a number of time steps in the 
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past. In inner product MGs, in contrast, the external information at stage i will consist of a 
p-dimensional vector (generalizing the 'information vector' of [14| |T5]), where a = p/N 
remains finite as N ^ oo, with entries 

h{e) = /[(I - C)A{i - A) + CZ{£, A)] (2) 

and with X £ {1, . . . ,p}. The function /[.] (which we choose to be anti-symmetric) allows the 
information given to the agents to involve more general properties of the overall bids than just 
their sign, on which the original model [HIS] was based. The Z(i, A) are zero-average Gaussian 
random variables, which represent the 'fake' alternative to the true market data A(£ — A), and 
p is the number of iteration steps in the past for which market information is made available. 
For reasons that will become clear we write = p^^J-\[l, A, Z], so 

j^x[e, A, Z] = /[(I - OM^ - A) + czii, A)] (3) 

If we ensure that /[.] remains finite, one has limTv^oo Jlxi^xi^)]'^ — ^'(1) for all i. Definition ([2]) 
implies that for C = 1 the external information is fully random, whereas for ^ = it represents 
true historic market data. Following pT] we distinguish between two types of definitions for 
the statistical properties of the 'fake memory' variables: 

consistent: Z{i, X) = Z{£ - X), {Z {i)Z {£'))= S'^6u' (4) 

inconsistent: Z(^, A) all independent, {Z{i, X)Z{i', X')) = S'^du'dw (5) 

We take the agents' strategies to be represented by vectors R*'' G IR^ with components In 
accordance with standard definitions |14l [T5] we choose these components to be independent 
zero-average and unit-variance Gaussian random variables, assigned before the start of the 
game, and remaining fixed throughout. They represent quenched disorder. The arrival of 
information at step i prompts each agent i to submit the bid bi{i) = ^\Rx'''^^^ Ix{t}^ 
where ai{t} denotes the preferred strategy at stage £ of agent i. In order to detect which of 
their S private strategies to use, all agents keep track of valuations Pia{i) of their strategies. 
These measure to what extent each strategy would have led to minority decisions if it had 
been used all the time: 

P^a{^ + 1) = Ka(^) " ^ A[l) ^ E^^hii) (6) 

The factor 77 > represents a learning rate. The active strategy Oj(^) of trader i at 
stage (. is chosen to be the one with the highest valuation pia{(.) at that moment, i.e. 
Oj(£) = arg maxa{j»ja(£)}. Our full equations now become, in more explicit form and with the 
notation ([3|): 

Pia{i + 1) = Ka(^) - ^ E -^aI^, a, Z]A{1) (7) 

=;7i^EE<'^'^-^A[AA^] (8) 



A 



Henceforth we restrict ourselves to the case 5 = 2, two strategies per agent so a G {1, 2}, and 
introduce the conventional definitions 

as well as = A^^^/^ This allows us to write R**^*^^) = a^*+sgn[(7i(-^)]^*. We then proceed 

to include decision noise in the familiar way via the substitution sgn[gj(£)] a[qi{l), Zi{l)\, in 
which the {zj(^)} are independent zero average random numbers, described by a symmetric 
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and unit-variance distribution P{z). The function (y[q,z\ is taken to be non-decreasing in 
q for any and parametrized by a parameter T > such that (T[g, z] G {—1,1}, with 
hmy^o <^\q^ a — sgii[(7] and Ivcht-^oo J dz P{z)a[q, z] = 0. The two main examples are additive 
and multiplicative noise: 

additive : a[q, z] = sgn[q + Tz] (9) 

multiplicative : z] = sgn[q] sgn[l +Tz] (10) 

T measures the degree of randomness in the agents' decision making, with T = bringing us 
back to ai{£) = arg maxa{pja(^)}, and with random strategy selection for T = oo. In on-line 
MG theories one will in practice only need the decision noise average 

a[q] = Jdz P{z)a[q,z] (11) 

Upon also adding the usual external perturbation fields {Oi{i)} to define response functions, 
our microscopic equations (j7|8p are then replaced by 

q,{e + 1) = q,{e) + 9,{i) - -JL= ^2 Cim, A z]A{e) (12) 
A{e) = A,{£) +^ {n^+^Y.''Me),z,{m\}m,Az] (13) 



The values of A{£) for i < and of the qi{0) play the role of initial conditions. 

Equations (|12ll3p are identical to those studied in [11] for look-up table MGs; the only 
difference between look-up table and inner product definitions for the information-to-decision 
conversion in the MG is in the definition of the function >1, Z]. We may thus take over 
from [11] the derivation of the effective single agent problem, and simply make the appropriate 
substitutions in the final resuW 



look-up table: Ag{-1,1}*^ 2^^^ = p, J'^^[i, A, Z] = VaN 5x^X(^e^A,z) (1^) 

X^{£,A,Z)=sgn[{l-C)A{e-^,) + CZ{£,fi)] 

inner product: Xe{l,...,p}, J'xil, A, Z] = /[{!-() A{e- \) + CZ{e, X)] (15) 

The inner product models of [14:\ [T5] correspond to choosing recipe (|15p . with C = 1 (fake 
history) and f[A] = A, giving Tx[i, A, Z] = Z{i, A). 



3. Generating functional analysis of generalized MGs 

3.1. The effective single agent equation 

We choose the real-time duration of individual iterations of the MG equations to be 5^ = fj/2p. 
This was proposed within the replica approach [5l[6l[T0], and later confirmed to be the canonical 
scaling via generating functional analysis studies [9l [H] . The generating functional analysis 
in [11], which via (fT5]) can be made to apply also to inner-product models, leads as always to 
closed deterministic dynamical order parameter equations which are fully exact in the limit 
N ^ oo and on times which do not scale with A^. The order parameters are the disorder- 
averaged single-site correlation and response function, which in terms of the original discrete 

I In fact, the formulation of the derivation in 'llj was chosen with the present generalization in mind. 
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iterations £ = 0,1,2, .. . are defined as 

C {£,£') = Jim (16) 

I 

G(A f ) = ji^n^ i E ^ (17) 

The brackets (...) refer to averaging over all realizations of both the decision noise variables 
{zi{l)} and the fake bid variables {Z(^, A)}, at all times. It was shown in [11] (to which we 
refer for details) that for N ^ oo and in terms of new time variable t = £5^ (which will become 
real- valued in this limit) the order parameters (|16|17p are to be extracted self-consistently from 
the following disorder-free effective single agent process with a retarded self-interaction and 
zero-average Gaussian noise rj{t) with non-trivial temporal correlations {r]{t)'i]{t')) = T,(t,t'): 

d /■* 

j{t) = e{t) -a dt' R{t, t') a[q{t')] + ^ r]{t) (18) 



dt 

As a result of the limit ^ cx) we must also write C(^,f ) C{t,t') and G{£,£') G{t,t'). 
One always has G{t, t') = for t < t' (due to causality), and C{t, t') = C{t', t) with C{t, t) = 1. 
If we write averages over realizations of the non-Mar kovian process (jlSp as (.••)*, the kernels 
{C, G} must satisfy for t > t': 

C{t,t') = {a[q{t)]a[q{t')]), G{t,t') = ^(afe(t)]). (19) 

The relation between the auxiliary kernels {R, S} and the order parameters {C, G} was found 
to be defined via an effective process for the overall bid A{£), namely by 



(20) 

i=t/5N/'=t' /5n 



nt,t')=fi \\^^ 5n-' {{w[£,£';{A,Z}]A{£)A{£',^^^^^^^ 



(21) 

l=t/5N/'=t'/SN 



Here ((. . .))a,z refers to an average over this effective stochastic process for the bids {^4} and 
the pseudo- history {Z}, a process that was found ^1] to take the form 

A{£) = Aei£) +cPe-hj2 Gi£, £'; {A, Z}]A{£') (22) 

with zero-average Gaussian random fields {</>}, characterized by 

{Mt){m,z} = \[l + C{£,l')]W[£,£'; {A, Z}] (23) 

The function .] in these formulae is defined as 

W[£,£';{A,Z}] = -Y^Ty\£,A,Z\TyX£',A,Z\ (24) 
Pa 

The macroscopic theory is closed. The remaining problem is to (i) express the kernels {i?, E} 
in terms of {C, G} using (j20l2ip . followed by (ii) calculating the order parameters {C, G} from 

(USD. 

It is quite satisfactory to find that the theory of [11] can be made to apply to a much 
broader class of MG models, which saves us from having to redo the generating functional 
analysis here. Mathematically, the differences between the two main MG families (look-up 
table versus inner product strategies) are within the macroscopic theory found to be limited 
to which expression to substitute in equations ()20|21|22p for the function W\£, £'; {A, Z}], which 
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measures the similarity between the market 'histories' (whether fake or real) as observed at 
stages i and i' of the process: 

look-up table : W[£, {A, Z}] = 5A(f,A,z),A(«',A,z) (25) 

inner product : W[l, i'; {A, Z}] = (26) 
^ ^ /[(I - OA{i - A) + CZ{i, A)] /[(I - OA{i'- A) + CZ{i', A)] 

X 

Here X{£, A, Z) E {—1, 1}^''^ (with 2^^ = p) denotes the 'history string' as observed by agents 
in the lookup table MGs at time i, with entries Xk{i, A, Z) = sgn[(l - C)A{i - k)] + CZ{i, k)]. 
Similarly, the differences between the two 'fake memory' definitions ()4|5p are limited to the 
details of the zero-average Gaussian variables {Z} in ()25|26p . and the associated averaging 
process (. . .)z- 



3.2. Time translation invariant stationary states 

In the special case of fully ergodic and time-translation invariant states (TTI) without 
anomalous response, where C{t,t') = C{t - t'), G(t,t') = G{t - t'), R{t,t') = R{t - t') and 
T,{t,t') = $](t — t'), one can derive from the effective single agent equation relatively simple 
and familiar expressions for persistent order parameters [TT]. The relevant scalar quantities 
are x = Jo°dt G{t), XR = Jo°dt R(t), c = limt_*oo C{t), and the fraction cj) of 'frozen' agents 
in the game. They were found [llj to obey: 

= 1- Erf M (27) 
c = [oo] { 1 - Erf [n] + -^Erf M - ^e""' } (28) 



2u'^ ' ' u^y^^ 

X = ErfM/axR, (29) 

with the short-hands u = ^/ax^(^[oo] /S0V2 and = S(oo). In order to find the TTI 
stationary solution {(/>,c, x} and the phase transition point (defined by x ^ cxd), we therefore 
do not need to solve for our order parameter kernels in full but just need to extract expressions 
for Xr and Sq from the stochastic overall bid process ([22|) . These latter two quantities can be 
written as 

wm{A,z}] +E^((^[^>o;K^}]^(^)»,^,,J (30) 

L 

si = lim hm {A, Z}]A{l)Aii'))) (31) 

Note that A(fi) and VF[0, 0; {A, Z}] express the initial conditions; for the purpose of evaluating 
the stationary state we may now drop the global bid perturbations {Ae{t)}. 



3. 3. Expressions for the kernels R and S 

In order to work out the theory also away from TTI stationary states and compare it to that of 
look-up table MGs, knowledge of just Xr and is no longer sufficient; one needs to calculate 
the kernels R and E in full. Here it is no longer clear to what extent the analysis of [llj can 
be adapted. The differences between the two model families start to manifest themselves. We 
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A,+i(4, ...,ir) = p'{(w[io, ir-, {A, Z}] n W[i,-i, li- {A, Z]] 



i=l 

r+r'+l//TXr\0„ B'-SA 7\lW 



{A,Z} 



[ n w[^,,u^^■. {A, z]]] [ n 4; {A, z}] 



{A,Z} 



(32) 



(33) 



As was done in in [llj we re-write the global bid equation (j22p as 

^ [dw + \vG{l,e)W% e- {A, Z]\]a{1') = Ae{t) + (^i 

and we invert the operator on the left-hand side, using 5n = ri/2p: 

A{t) = A,{i) + <l)t + Y,{-5NPY G{£,ei)G{£i,e2)...G{£r-i,er) 
r>o ei...er 

X W[e, h\ {A. z}W\hA2\ {A. z}]... W[ir-l,ir; {A, Z}] [AM + </>^] 

We insert (j34p into (|20p . and consider only infinitesimal external bid perturbations A^: 
R{t,t') = 6{t-t')+ lira ^\y2{-6NY V G(4,£i)G(£i,^2) • • • G(^r-i,^r) 



(34) 



P^Y(h^[£o, 4.; {A, Z}] n {A, Z}] 



i=l 



{A,Z} 



1 



,5(t-0+ lim — <^5](-<W)'^ E G(£o,^i)...G(£,_i, 

(V^O Oat ' 



r>0 



X Ar+i(4, . . . ,^r) 



^0 = 



(35) 



This is clearly identical to the corresponding expression found in [llj (although now the 
meaning of the kernels Afc(. . .) is allowed to be different, dependent on which model family we 
choose to apply the generalized theory to). Similarly we can insert (|34p into (j2ip . again with 
Ag — > 0, and find 

S(t,0=r/ lim -l| {-^nY^''' E G{£o,ii)...Gilr-i,lr) 

X J2 G(^[„^'i)...G«,„i,<0 {^eJe',)i^\A,z}W[io,i'o-AA,Z}] 



[ n W[t,.^,t,; {A, Z}]] [ n W[e'^^„t'^; {A, Z}] 



i=l 



{A,Z} 



hm <^ Y (-^nY^"' E G{£o,ei)...G{£r-i,ir) 

X J2 G{£',,e[)...G{e',,_„e[.)[i + c{£r,i'r>)] 



(36) 
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The limits (5/v — > in (j35l36|) are well-defined. Each time summation combines with a factor 
5i\f to generate an integral, whereas pairwise identical times in (j36p leave a 'bare' factor 5^ but 
can also be anticipated to cause Aj.+r'+2{- ■ •) gaining a factor p = fj/25]\f in compensation. 

Equations (j35|36p show, upon inserting the respective definitions (|25l26p of the function 
VF[.;.], that replacing lookup table strategy definitions by inner product ones does have 
implications. For instance, the following relation holds for lookup table MGs but appears 
to be no longer valid for inner product ones: 

Ar+r'+2(^0) ■ ■ ■ , ir', i'o, ■ ■ ■ , ^'r') = ^r+r'+2(^0) ■ ■ ■ , ir , ^0, ■ ■ ■ , ^'r') (37) 



4. Inner product MGs with inconsistent fake market history 

The simplest instance of our inner product MG model is the case C = 1) i-e- fake market history, 
of the inconsistent type. Here the function W[iJ';{A,Z}] = p-^ Y.x f[Zii, \)]f[Z{i' , X)] is 
no longer dependent on the real market bids A, and {Z{i, X)Z{i' , X')) = S'^5ui5\\i. Even for 
this case no exact macroscopic solution has so far been published. We define the short-hands 

~ 112 

Zt^x = f[Z{l,X)\ and k„ = ^Dz /"[Sz] with the Gaussian measure Dz = (27r)~2e~2^ dz 
(since f[A\ is anti-symmetric, = for all odd n). 



^.i. The retarded self-interaction kernel R 

To proceed we have to calculate the two functions ()32|33p . They occur only in expressions (f35]) 
and (I36p . where causality of the response function G enforces a helpful ordering of the time 
arguments. The kernel Ar+i(^0) . . . ,^r) as occurring in ([35]) is easy to evaluate, since here we 
may use £o > > • • • > ^r- This property, in combination with ki = {Z(_^x) = for all (£, A), 
implies that the only non-zero contributions are those where the Gaussian variables {Z} have 
pair- wise identical indices: 

Ar+l(4, ■■■■>^r) = - E {Zi^^XoZirM n -^^«-i,A,^^>,A,) 

" Xo...Xr = l i = l ^ ' 



E {Z£oMZeoM)[ Y[{Ze„x,Zi^^x,+i)){Ze,.,XrZe,,Xo) 



V 

^ Ao...Ar=l i=l 

P r-l 



K 



r+l 
2 



^ Ao...A,,=l i=l 



E KM{\{K^.+.)^^rM=^l^' (38) 



Insertion into (j35p then gives us the fully explicit form 

R{t,t') = K2{5{t-t')+ \\m ^Y^^-SnYkI V G{to,h)...G{lr-l,lr)\ 

= K2[5{t-t')+Y{-K2YG'{tX)} 

r>0 

= i^2{^ + K2G)-^{t,t') (39) 

This result for inner product MGs depends on the function /[.] and the variance S only via a 
single parameter K2 = / Dz /^[S'z], and reduces to the corresponding expression found earlier 
[9] for look-up table MGs with fake histories when K2 = 1. 
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We now turn to the function (j33p . This is needed in ()36p but only for time combinations with 
io > £i > . . . > £r and i'^ >£[>...> i'^,. Again we first group together the various terms 
that have identical time labels, and we also introduce two auxiliary summation indices A that 
allow us to separate the terms with times of the type {£{} from those of the type {i'j} in a 
clean way: 



^r+r'+2(4, . . . ,ir;£o, ■ ■ ■ ,4') 
P P 

E E 

Ao...A,.=l A'...A',^ =1 

1 r' + l 



~ E E (^i;o,Ao^^;,,A„%,A',^,^^',,A', , 



[n4-.,A.M[n%_„A^^. 
i=\ j=i 

^ E E '^A,+ i,A;,^/A(„Ao(^4,Ao%,Ai^f.,A.+ i%„A;,,^^ 

r r' 



^ Ao...a,+i=ia;,...a;,^^=i 



^£'.,A'. 

=1 i=i 



{^} 



E E '^A.+ i,A;,,^/A;,,Ao([n4A4,A.+i][n%,A;%,A;+J)^^j (40) 



^ Ao...A,+i=i a^...a;,^^=i ' 'i=0 'j=o 



Due to the time ordering relations of the problem, there can never be time coincidences 
amongst members of the set {io, . . . ,ir}, nor amongst members of the set {^g, . . . , Thus 
the various factors within the first pair of round brackets ( 111=0 • • • ) in (|40p are all independent, 
and so are those within the second pair of round brackets ( 115=0 • • • ) • only possible time 
coincidences are pairings between un-primed times {£+ and primed times {i'j}- Thus we have 
to consider the following cases: 

• If there are no time coincidences, i.e. ii ^ i'j for all the evaluation of (|40p proceeds 

similarly to that of the previous functions Ar+i{. . .): 

Ar+r'+2(^0; • • • ;-^r;^0, • • • , C') 

P P 



_ 1 

— — ^2 

P Ao...A,+i=i aj,...a;,,^^=i i=0 j=0 

= I -r'^' E E '^Ai.Ao = -r'^' (41) 

f Ao=i a;,=i 

Now consider the effect of a time pairing, where = I'j (note: their can be multiple 
pairings, but the number of coinciding times is two at most, due to the built-in time 
ordering). The contribution of the Gaussian averages with times to the sum 

J2xo. Xr+i 5I]a(, ..a', (00])) which in the absence of pairing equaled simply 

no pairing : (Z^^,a,^4,a,+i^£;,a^^^^.,a;^,) = (^4,A,^f.,A,+i)(^f;,A^^£;,A^^,) 
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no pairings: <5a„a,+i5a'.,a'.^^ 



j J + 1 



r' + 1 



h„x,+^Sx'^,xr^^ ^a„a;^/a,+i,a; 




j + 1 



l'^A' A 



J + 1 



r' + 1 



1 rr+l 

A = ^j-- ^a.,a.+i^a;,a;^_^^ ^ Sx^^x'A,,x^+i^x'.,x'^^^ 

1 

Figure 1. Diagrammatical representation of the different contributions of the 
Gaussian averages to the function Ar+r'+2- Each label i S {0, . . . ,r + 1} and each 
label j G {0, . . . , r' + 1} is drawn as a distinct vertex of a graph. Each factor Sx^^y. 
is drawn as a line segment connecting the vertices i and j. Top graph: the case 
where there are no time coincidences. Bottom three graphs: the three different new 
contributions that are generated by the occurrence of a time pairing where £i — i'^ . 




now becomes: 

ii = £'j : (Z^^,a,^4,a.+i^£;,a;%,a;_^^> = (Ki - 3kI) 5a„a,+i5a;,a;^/a„a; 

+ 4 [^x,,x,+i5xr,xr^, + ^\,y,^\+i,y,+, + '^a„a;+/a.+i,a;] 

The effect of time pairings can be represented diagrammatically; see figure [H 

We now make the crucial observation that in expression (|36p any time pairing would 
inevitably generate a 'bare' factor 6n that would no longer be absorbed into an integration 
via X^^fc — > Jdtk- Hence the only contributions to Ar+r'+2 to survive the limit 5]y ^ 
in (j36p are those where each time pairing also generates an extra 0{p) factor to compensate 
for the emerging 6^- Let us inspect the various diagrams, see figure [T] (and the higher order 
versions generated due to multiple time coincidences) and their contributions to Aj.+r'+2- 
Each connected diagram implies that ultimately in the summation X^Ao.-.Ar Saq...a', 
restricted to Aq = Ai = ...A,. = Aq = ... = A',/. This leaves a sum over Aq only, and a 
final 0{1) contribution to Ar+r'+2- On the other hand, in the case of a disconnected diagram, 
each internally connected sub-diagram will give an 0{p) factor. Here the order of the final 
contribution to Ar+r'+2 is where L denotes the total number of time pairings that gave 
rise to the diagram cuts. Hence the only relevant diagrams in figured] are the top one (when 
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there are no pairings) and the third from the top (when there are pairings, but where these 
also generate compensating 0{p) factors), including the higher order diagrams with multiple 
vertical connections (in the case of multiple time pairings). 

It follows, in combination with the contribution coming from the unpaired terms and 
from each possible time pairing as calculated earlier, and taking into account the crucial time 
orderings in (136]) . that 



i=lj=l 



r+r'+2 



nn 

1=1 j=i 



1 + ^6,^^,,{1 + 0{5n) 



(42) 



Insertion into (j36p gives 

^{to, t'o) =4Y. (-'^2)"+" dti... dtrdt[ . . . dt',, n n 



r,r'>0 



1 + l^V^iU - t'. 



(43) 



i=lj=l 

X [1 + c{tr, t;,)]G'(to, ti) . • . G(t,_i, t,)G(t'o, t'l) . . . G(t;,_i, t;,) 

We see again, as with the retarded self-interaction kernel ii, that also this result for inner 
product MGs reduces to the corresponding expression found earlier [9] for look-up table MGs 
with fake histories when K2 = ^■ 



4-3. Summary, TTI stationary state and phase diagram 



The full and closed dynamical equations for the inner product MGs with inconsistent fake 
market information are thus found to be given by 

6 



C{t,t') = {a[q{t)]a[q{t')]). 



G{t,t') 



66 {t 



(44) 



Averages are defined with respect to the effective single agent process with a retarded self- 
interaction and an effective zero-average Gaussian noise r/(t): 

d 



dt 



9(t) - aK2 

Jo 



(45) 



ivitoMt'o)) = 4J2 i-^2Y^" / • • • dtrdt[ . . . dt',, n n 1 + ~ ^'i^ 

r,r'>0 ° 4=1 i=l 

X [1 + C{tr, t;,)]G'(to, ti) • ■ ■ G(t,_i, t,)G(t'o, t'l) . . . G(t;,_i, t;,) (46) 

We notice that the role of the remaining parameter K2 is to define a characteristic time scale 
for the process, and that apart from this our equations reduce in the limit of infinitesimal 
perturbation fields exactly to the theory of the fake history look-up table MGs as derived in 
[9]. To see this we re-define G{t,t') = K2^G{t,t'), 9{t) = K29{t) and 7]{t) = K2r](t), and find 

= {a[q{tMq{t')]U G{t,t') = ^(a[g(t)]). 



C{t,t') 
Id,, 

-^2dt''^'^ 

mo)m'o)) 



e{t') 



§{t) - a Cdt' (I + G)-^{t, t') a[q{t')] + ^r)(t) 
Jo 

{-ly^'' / dti... dtrdt'^ . . . dt',, n n 

r,r'>0 ° i=l 7=1 



(47) 
(48) 



1 + :;VS{ti - t',- 



X [1 + c{tr, t',,)]G{to,ti) . . . G(t,_i, t,)G(t[), t'l) . . . G(t;._i, t;o 



(49) 
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It will no longer come as a surprise that upon using ([c 



to calculate in time-translation 
invariant stationary states the persistent order parameters = Jdt R{t) and Sq = S(oo), we 
find that the only difference between the stationary state order parameter equations of lookup 
table MGs and the present inner product MGs is a re-scaling of the static susceptibility x- 



u = a[oo\\J a/2{l + c) 
= 1- Erf [u] 

c=a\oo]{l- ErfM + -^Erf^ - ^e-«'| 
X = Erf[n]/K2(a-ErfM) 



(50) 
(51) 

(52) 
(53) 

Since phase transitions in the present type of MGs are defined by a divergence of x, we conclude 
that not only the values of the static observables {0, c} but also the phase diagrams of the two 
fake history MG model families are identical. 

4.4- The volatility 

Finally we calculate the volatility. Since the average overall bid is zero, the volatility is here 
defined by cr^ = limt^oo hm5jy^o((^^(^))){A,z}l£=t/5jv • Inserting into the definition of cr, 
with Ae 0, gives 

a^= lim hm \ ^ (-5^)'^+^' ^ G{io, h) ■ ■ ■ Giir-i, Q 

X G(£^),ri)...G(<,,_i,£;0'^^A'5A^^P'''''^'(( ^'^^'■'^C'^MA^} 

X [ n W[ii.u^u {A, Z]]] [ n Tr[4_i, 4; {A, Z]] 
i=i j=i 

= hm hm I I Y (-^nY^^' E G(4 



{A,Z} 
h)...G{ir-l,er) 



X ^r+r'+li^O, ■ ■ ■ ,^r',^0, ■ ■ ■ 



with 



Aj,_(_r'+l(^0) . . . ,-^r;^0, • • • ,£'r') 

= p^+^'{{W[£r,KrAA, Z}] [ n {A, Z}]] [ f[ W[£'^_„ey,{A, Z}] 



(54) 
(55) 



{A,Z} 



i=l j=l 

We need this function (fSSl) only for times with Iq > £i > . . . > £r and £q > £[ > ... > £'^,. 
For the case of inner product MG, where W[£,£'; {A, Z}] = p^^ Sa=i Ze^\Zi'^x, its evaluation 
is very similar to that of ([33|) : 

_ p p 

^r+r' +l{io, ■■■ Jr;£o, ■■■ ,£r') = - E E "^A^+l.A', 

^ Ai...Ar.+i=l A'...A',^ =1 

i r' + l 



^^o,ai^^[„a; [n ^«nA,^4,A.+i] [ n %-a;^. 



£'.,A' , 



{^} 



(56) 
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Due to the time ordering, the factor Zi^^x^Z^/^ y^ is statistically independent of all others, so 
that 

1^2 



Ar+r'+l(^0, • • • ,-^r;^0, • • • ,^r') = X! X! 

^ Ai...A,.+ i=l A'...A', =1 

1 r' + l 



'^A.+i,a;,^^ 



7 

[n 

1=1 



Zf, X, Ze 



Again we inspect the effect of time pairings: 

• If there are no time coincidences, i.e. ii ^ i'j for all we simply find 

\-+r'+l{^0, ■ ■ ■ Ar'A'o, ■ ■ ■ A'r') 
1 



r+r'+l 



P 



P P 

E E 

Ai...Ar+i=l A'...A',^ =1 

1 r' -\-[ 



'^A,.+i,A' jn'^A„A,+i][n'^. 
i=i j=i 



(57) 



^.+.'+1 ^ ^ 



Ai 



K 



r+r'+l 



(58) 



Ai=l A'=l 



The effect of a time pairing ii = i'j on the contribution of the Gaussian averages with 
times is exactly as before: 

no pairing : (Z^^^a.^^^a.+i^^^.a^^.a^^,) = (^^„a,^£„a,+i)(^^;,a;^^^,a;^,) 



J' i+i 



Ii = e'j : (Z^,,A,^4,Ai+i^£^,,A;^^;,A;_^j> = («^4 - 3^2) '5a,;,A.+ i'^a;,a;^/a„a; 

+ ^2 ((^a^a.+i^a^.a;^^ + ^a^a^-^a.+i.a;^^ + 5a„a;^/a,+i,a;) 

The corresponding diagrammatic representation is shown in figure [2j Since each time pairing 
will have to be compensated by an 0{p) factor to retain relevance in the limit 6]\[ 0, the 
only relevant diagrams in figure [2] are again the top one (no pairings) and the third from the 
top (a time pairing with vertical connections that cause a diagram cut), including the higher 
order diagrams with multiple vertical connections (in the case of multiple time pairings). We 
conclude that 

r r' 

A,,+,„+i (4, ...jr-jo,---, I'r') = 4+'''+' n n [1 + h/', (p + ^ip") 

1=1 j=l 



^r+r'+l 



nn 

i=ij=i 



(59) 



and hence 



cj2 = -K2 lim y (-K2Y~^'' / dti... dtrdt[ . . . dt', TT n 

r,r'>0 ^ 1=1 j=l 



X [1 + C{tr,t',,)]G{t, ti)G{tut2) . . . G{tr^utr)G{t, t[)G{t[,t'2) . . . G{t',,_M 



^lirn — Cdt J2 (-^2)'+"' Hdti . . . dtrdt[ . . . dt',, f[ f[ 
-r^oo 2r Jo ry>o ° i=ij=i 

X G{t, ti)G(ti, is) . . . Gitr-i,tr)G{t, t[)G{t[,t'2) . . . G(t;,_i, t;,) 

X [1 + C{tr,t',,)] 



1 + l^¥iti - <•) 



(60) 
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no pairings: <5a„a,+i5a'.,a'.^^ 



J' j+i 

1 2 j j + 1 r' r' + 1 



h„x,+^Sx'^,xr^^ ^a„a;^/a,+i,a; 




j + 1 



l'^A' A 



'^a„a;'5ah 



J + 1 



r' + 1 



1 2 i i + 1 r r + 1 

A = ^j-- ^a.,a.+i^a;,a;^_^^ ^ (^a^a^^a^a.+i^a^.a;^^ 

1 2 

Figure 2. Diagrammatical representation of the different contributions of the 
Gaussian averages to the function A^+r'+i- Each label i S {1, . . . ,r + 1} and each 
label j e {1, . . . , r' + 1} is drawn as a distinct vertex of a graph. Each factor ^a^.a' 
is drawn as a line segment connecting the vertices i and j. Top graph: the case 
where there are no time coincidences. Bottom three graphs: the three different new 
contributions that are generated by the occurrence of a time pairing where ii = £'j . 




Comparison with the corresponding expression in [9j for look-up table MGs shows that the only 
difference is in the various occurrences of K2 ■ If we implement the transformation that maps the 



effective single agent exactly onto that of the lookup-table models, viz. G{t, t') 



'G{t,t'), 



we find that the inner-product volatility differs from that of the lookup table models by a 
factor K2: 



a 



lim — 



dt 



E 

r,r'>0 



(-1) 



r+r' 



dti 



dtrdt[ . . . dt'^, W W 
i=ij=i 

X G{t, ti)G{ti,t2) . . . G{tr-utr)G{t, t[)G{t[,t'2) . . . 

X [i + c7(t,,,t;,)] 



1 + l^iti 



(61) 



5. Inner product MGs with real market history 

The previous cases could be solved exactly and in full, due to the absence of real history. 
We now turn to the more demanding situation where ( < 1 in definition (j26p . so that the 
information vector truly depends on the past market history. Here it will turn out advantageous 
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to define the stochastic and time dependent symmetric p x p matrices B(£) with entries 

Bxx' (i) = Tx A, Z]Fx' [I, A, Z] (62) 
with ^;^[. ..] as defined in ([3]). This allows us to write the relevant functions Ar+i(...), 

Af.+r-'+2(- • •) a-iid A,._|_r/_|_i(. . .) (that occur in the kernels R and S and in the volatility) as 
averages over a trace: 

A,+i(4, . . . ,4) = ^Tv((b(4)B(£i) . . . B(4)))^^^^^ (63) 
A,+,,+2(4, . . . , . . . , Q = ^Tr((B(4)B(£i) . . . Wr)Wr') ■ ■ ■ B(^i)B(^o)))^^_^^ 



Kr+r'+i{(oJi ■ ■ -Jr-JoJ'i, ■ ■ • ,C) = -Tr ((b(4)B(4) • • • B(4)B«,) . . .B(/i 



These expressions are of course still exact, but not easy to evaluate. 



(64) 
(65) 



5.1. Short history correlation times in TTI stationary states 

We now make an approximation for the TTI stationary state solution in the spirit of the 
short history correlation times ansatz made in [llj for lookup table MGs; in fact we will find 
below that for lookup table MGs the two are identical. In working out our expressions for 
XR = Jdt R{t), S(oo) = limT-^oo ^{t + T,t) and the volatility o"^ we will replace the kernels 
(j63|64l65p by the values they will take for times which are sufficiently widely separated to 
de-correlate the random matrices B(£), i.e. we replace 

f[B{ei) ^ B'{A,Z) B{A,Z)= Ihn ]rJ2B{£) (66) 

^=l ^ I^L 

For TTI states this allows us to express everything in terms of the average eigenvalue 
distribution g(fj,) of the now time independent random p x p matrix B(A, Z): 

Ar+r'+2(^0i • • • ,^t',^0: ■ ■ ■ 



Ar+r'+l(^0;-^l ■ ■ ■ , ^r', ^0, ^l, ■ ■ ■ , ^'r 



Since B{A, Z) is non-negative definite one always has /x > 0. Insertion of (|67l68l69p into 
expressions ()35|36p for the kernels R and S followed by appropriate integration and re- 
summation of the series leads us to the following expressions for xr and Sq = S(oo): 

POO POO 

XR= dfi Q{fi)— Si = {l + c) dfi Q{fi)— "2 (70) 

Jo 1 + /^X ^0 (1 + I^X) 



= -Tr ((B'-+i(A^))){A,z} 

POO 

= / dfi Qif,) (67) 
Jo 

= ^Tr{{W+^'+\A,Z))){A,z} 

POO 

= / df, Qif,) +^ (68) 
Jo 

,) = lTr ((B'-+^'+1(A^)))m,z} 

POO 

= / d^^>(M)^^+'■+l (69) 

JO 
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The effects of having real memory are concentrated fully in the eigenvalue distribution g{fj,). 
Once this distribution has been calculated in terms of {c, (p, x}, the persistent order parameters 
of the MG will be given by a closed set of equations. Upon introducing a convenient parameter 
a; G [0, 1] via 

Jdfi g(^)/XX/(l + /iX) fry,. 

^Jdfi Q{fj,)[nx/{'^ + fJ'X)? 
our closed order parameter equations take the form 

.= ""[°"]^, W^r-^^^^ <^ = l-ErfM (72) 
v/2(l + c)' a Jo l + (/Ux)-i' ^ ^ ^ ^ ^ 

c = [oo] 1 1 - Erf [d + At Erf [n] ^e""" | (73) 



From this we can already extract several results, and explain observations made in the past 
for MG types other than look-up table ones on the basis of numerical simulations: 

• Irrespective of the extent to which the histories are real, and irrespective of the nature of 
the function f[A], the ergodicity-breaking phase transition point marked by a divergence 
of X occurs always at the value adT) that was found originally for fake history look-up 
table MG models [3 El [71 IB] > as solved from the closed set 



u = o-[oo]^a/2(l + c), Erf[n] = a, (74) 
c=a^[oo]{l- Erf [n] + -^Erf M - ^e~"' | (75) 



• Equations (|71|72|73p return to those derived earlier for inconsistent random histories 
when uj = 1, whereas real histories affect the persistent observables as soon as w < 1. For 
finite x we thus arrive at the conclusion, via ()7ip . that our observables are affected by 
having real histories if and only if g^fi) is not a ^-distribution. A nonzero width of g{fj,) 
is the fingerprint of real histories influencing the behaviour of the MG. 

We note that the combined equations ()70|71|72|73p are identical to those found for the lookup 
table MGs with true market histories [llj, with the eigenvalue distribution q(^) taking over 
the role of the relative history frequency distribution in [llj. This can be understood. 
Our choice of notation allows us to apply the present derivation also to the lookup table 
models, simply by making in (j62p the substitutions A G {!,..., p} — > A G { — 1,1}*^ 
and J^x[i,A,Z] —>■ ^ ^W{t Azy '^here X{i,A,Z) is the history string (potentially partly 
fake) as observed by the agents at step i. In lookup table MGs one would therefore have 

Here it^{A, Z) denotes the stationary state frequency of occurrence of history string A. Since 
this matrix B(A, Z) is of a diagonal form, its eigenvalues are simply the diagonal entries 
pTT^{A, Z), and the asymptotic eigenvalue distribution becomes 

gifi) = lim - Y,{{5 [/i - pTT^iA, Z)])){^^z} (77) 
A 

Thus, for lookup table MGs the eigenvalue distribution g{fi) is indeed identical to the relative 
history frequency distribution of [llj . 
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The question of whether and to what extent having real as opposed to fake market 
histories affects the MG thus boils down to assessing the shape of the eigenvalue distribution 
g{n). For look-up table MGs this distribution is known to become nontrivial above the critical 
point \10\ lllj. Numerical simulations show that the same is true for inner product MGs, so 
also here models with real history will behave differently from those with fake histories. 



5.2. Calculation of eigenvalue spectrum g{fi) for weakly correlated bids 

In the remainder of this paper we will deal exclusively with the inner product version of the 
MG with real history and inconsistent bid noise, i.e. from now on we will have 

Bxy{A,Z) = jJ2 f[{l-C)Aii-X)+CZ{£,X)]f[il-C)Aii-X')+CZ{i,X')] (78) 

where {Z{i, X)Z{i' , X')) = S'^6ui6x\', and where we focus on the regime where 1 <C p <C L. 
We try to understand the origin of the nontrivial eigenvalue distribution Q{fj,) in inner product 
MGs and its impact on the order parameter equations. Numerical simulations show that, even 
when the histories are real, the bid correlations are still very weak, see e.g. figure El this we 
will exploit in our analysis. The remaining complications arise from the fact that, although the 
matrix B(A, Z) is an average over a diverging number of at most weakly correlated variables, 
the number of entries of B(^, Z) also diverges. To prepare the stage we first write the time- 
translation invariant bid covariances as H(A) = {A{£)A{i + A))|^^2}., and introduce the short- 
hands 

Qo= Jdz /2[^y'(l_02^2 + ^252] (79) 
Qi = fDy\ [dz f[{l-(:)ay + CSz]Y (80) 



with = limg^oo{{A'^{^))){A,z} = -(0), and < Qi < Qo- As for the look-up table MGs, 
one expect also for inner product MGs the distribution Q{p) to be well-defined as p oo, 
since Jdfi Q{fi)fi = p'^ T,x<p{{Bxx{A, Z)))a,z = Qo- 

In the absence of bid correlations, as for fully fake histories, B(A, Z) will be self- 
averaging for L ^ oo and fixed p; this implies (due to the time-translation invariance of 
the bid correlations) that B will be of the Toeplitz form. One expects time translation 
invariance to hold also if the correlations between bids at different times are sufficiently 
small, and also B(j4, Z) to remain self-averaging at least with respect to the inconsistent noise 
variables, although possibly not with respect to the bids. Here we explore the consequences 
for the spectrum q(h) of assuming the bid correlations H(A) for A 7^ to be small, and 
of taking 'B(A,Z) to be Toeplitz and self-averaging with respect to the {Z{i,X)}, i.e. 
Bxx'{A,Z) = B{X — A'|yl). If H(A) decays to zero or to small random values sufficiently 
fast, we may continue B(74) L-periodically, so that it becomes a circular Toeplitz matrix. Its 
L eigenvalues //^ (r = 0, . . . , -L — 1) will now be given by 

L-l 

fir= J2 i?(A|A)e-2-^'^/^ = Qo + SriA) (81) 

£'=0 

with Qq as given in ([79l) and, upon replacing J2e<L by '12£<l dZS]) (which is allowed since 
L — > 00): 

Sr{A) = jJ2 e-'"'''^/^[/z?^ f[{l-OA{e-e') + CSz]] [Jdz f[{l-C)A{£)+CSz\ 



Minority games with inner product strategy definitions 



18 



f[A]=A: 



f[A] = sgn(A) 



C=o 



2000 -2000 





C=i . 





2000 -2000 



2000 -2000 






C=i 





000 2000 -2000 



Figure 3. Time-translation invariant bid covariances S(^) = {A{t + £)A{t)) of inner 
product MGs ([75)1 measured in numerical simulations, with N = 1025 and a = 2. 
Covariances are calculated over 200. time steps, following equilibration. Top row: 
f[A] A, with C e {0, i, 1} Bottom row: f[A] = sgn(A), with ( e {0, i, 1}. The left 
pictures refer to strictly real histories, the right pictures to strictly fake histories, with 
the middle ones representing a mixing of the two. These data are representative of the 
typical behaviour throughout the ergodic regime. The bid correlations are always seen 
to be weak, which justifies our approximate calculation of the spectrum g{fi). We also 
observe qualitative differences between f[A] = A and f[A] = sgn{A), e.g. vanishing 
volatility as C — > (strictly real histories) for f[A] = A, which will be addressed later. 



= 4f E /^^ /[(l-C)^W + C5z]e-2-V^ 

X 4f E /^^ m-OA{£')+CSz]e'-'''^/'^ (82) 

The first and second line are both sums over a large number of nearly independent random 
variables. If the correlations between the bids at different times are sufficiently small, we may 
take the second sum to be independent of i, which allows us to put = L in this term, provided 
we correct for the inappropriate generation in doing so of terms where the bids involve identical 
times (i.e. where i = £', as such terms were absent initially): 

SriA) = \ZriA)\^ - Ql (83) 

^riA) =^J1Jd' m-C)A{e) + CSz]e'^''^/'^ (84) 

Clearly ZL~r = Zr- Given our assumption of weakly correlated bids and given that f[A] is by 
definition anti-symmetric, the L complex random variables Zj. (where r = . . . L — 1) will for 
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L — > CO acquire zero-average Gaussian statistics (central limit theorem), with 

{ZriA)Zr, {A)) = ^ |j ,^^^Kr+r'yL ^ ^odL " (85) 

{MA)ZAA)) = ^Y.--"'''''^''^'^" = T^Sr-,r',mo,L " (86) 
{ZriA)zAA)) = ^J2^2.Uir-r'yL ^ ^ ^^^^ _ ^ (gy) 

If the spectrum is well-defined as L oo, we may take L to be odd and write the 
spectrum as an average over {zi, . . . ,-Zl/2-i/2}) using z^-r = 'Zr to deal with the remaining 
{•2L/2-1/2; • • • ) -^L-i}- We may disregard r = 0, since it gives only a vanishing contribution to 
the spectrum. The remaining -^^^ complex Gaussian variables are written in terms of their 
real and imaginary parts as Zr = Xr + it/r , for which one finds 

{XrXr') = 2(2^-1^^ ^"'' ~ {VrVr') = 2^1-1) ^'''''' i^rVr') = (88) 

If for L ^ cxD we are allowed to neglect the 0{l/L) correlations between the Gaussian Xr, 
we may use the law of large numbers and write the spectrum Q{fi) as an average over their 
distribution, giving (with \zr\'^ = x^ + y^)- 

^ 2 I „,2^ 



Q{fjL) = j DuDv 6[iJL -Qq + Qi- ^Qi('U^ + V 



It is encouraging to see that this simple expression passes the two tests at our disposal. Firstly, 
it obeys the exact equation / f-QifJ-) = Qo- Second, since lim^^i Qi = one always has 
lim^^i q{ij,) = 5[fi — Qo], so in the fake history limit it leads us correctly to the standard 
stationary state equations of the fake history MGs. The remaining deviations very likely to be 
finite size effects. Having a finite value of L in measurements removes the self- averaging of the 
matrix Bxx'{A, Z) with respect to the noise variables Z{i,X), which our L ^ oo calculation 
assumed (and thereby deforms the matrix in the direction of the random matrices as studied 
in e.g. \17\ \18 \ [T9]). and furthermore we have in these simulation experiments only been able 
to sample in equilibrium states matrices of limited dimension p € {1025,2049}. Both these 
limitations result in a smoothening of the spectrum as measured in practice, compared to 
the asymptotic shape. The effect of retaining the 0{L~^) correlations between the Xr in the 
above calculation turns out to be an order reduction of the width of ^?(/i); furthermore, 
when including 0{L~^) terms finds that the condition Jdjj, ng{n) = Qo is violated at order 
L~^, which implies (as expected) that in that order we can no longer treat the Zr as Gaussian 
variables. 

We note that the above reasoning would also apply for S(A) = a'^6xo, a situation 
indistinguishable from replacing the actual bids in the history signal by consistent fake 
information, so we must conclude that the origin of nontrivial eigenvalue distributions g{fi) 
is the consistency of the histories (although of course the precise shape of g{fj.) will depend 
on whether or not the histories are real). We note also that allowing for 'B{A, Z) not being 
self-averaging with respect to the bids has been crucial in our argument. Had we set the bid 
correlations to zero and averaged B(A, Z) over the bid process, we would have found the trivial 
spectrum g{fj,) = 5[fj, — Qo]. 



Let inspect what our general results imply for the two canonical choices for f[A]. For 
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H H fi 

Figure 4. Comparison between eigenvalue spectra g{ii) of the stochastic matrix 
B(A, Z) of inner product MGs ([75)1 as measured in numerical simulations (histograms, 
for N = 1025) versus the approximate prediction ([5^ based on assuming weak bid 
correlations, for f[A] ~ sgn(A). Top row: a — 1, with ( E {0,^,1}. Bottom row: 
a = 2, with C e {0,i,l}. All matrix entries in simulations were calculated over 
L = 400. time steps, following equilibration. For ^ = | we used in ((92)l the value for 
a as measured in the simulation; for C G {0, 1} the volatility drops out of the formula. 
Note that there are at least two obvious sources of deviation between prediction and 
experiment (apart from the crude approximation used in the theory): finite size effects 
{N = 1025 rather than N — > oo) and finite observation time effects {L < 400. iV rather 
than L — > oo). Given these limitations, the qualitative agreement between prediction 
and experiment is surprisingly acceptable. Spectra of random matrices are indeed 
known to have finite size effects which decay much slower than N~^/'^, however, we 
will find later that at the level of the MG's conventional order parameters the finite 
size deviations between theory and simulations are of the expected order N^^/"^. 





f[A] = A one obtains Qq = {I - Cfa^ + C'^S'^ and Qi = (1 - Cfa"^, so that 



(90) 



(i-C)V2 

Fully real history would give Vhuc^^q q{p) = 9[fi]a~^e~^^'^^ . For f[A] = sgn(^), in contrast, 
one has Qq = 1 and 



Qi 



Dy Erf' 



(l-C)tTj/ 

CSV2 



vr 



-arctan[W 1 + 2 



(91) 
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and so 

/[larctan[yi^^ii=<^]-l 

iarctan[^l+2ii=ig^]-l ^ ^^2) 

Here the limit of fully real history gives lim^^oQi = 1 and lim^^o = 0[/^]e~^. The 
persistent order parameter equations (|71|72|73p do not involve the volatility, but it turns out 
that we generally need to know the volatility in order to use (j92p and (|90p for closure, unless 
C G {0,1}, i.e. unless the histories are strictly real or strictly fake. In figured] we compare 
the above predictions with the eigenvalue spectra as measured in numerical simulations of the 
inner product MG with f[A] = sgn(^), where (for {0, 1}) we inserted in our formulae the 
value of a as measured. Given the severe finite size and finite time effects in such simulations, 
the agreement is surprisingly acceptable. 



;U-2+^arctan[A/l + 2 



2„2. 



—^TsT- 



/i-2+farctan[^/ 1+2^^=1^ 



5.3. The volatility 

We turn to the calculation of the volatility, which measures the fluctuations in the bids A{i). 
It turns out that these bids are always Gaussian random variables, see figure [SJ which is not 
obvious since the bids obey a generally nonlinear stochastic equation. The first step is to 
substitute ([65]) (obtained upon assuming short history correlation times in ergodic stationary 
states) into the as yet exact formula ([M]) for the volatility, giving 

t-ooV-oL2^,^g Jo ^^^^ 

X J2 G(^^„/i)...G(£;,_i,£;,) [i + ^^(^r,^;o]'^^A'^M,}L=iM ^^^^ 

t.g...t^, 

As always we next have to approximate non-persistent dynamical order parameters by 
expressions involving only persistent ones. Following procedures similar to those used in the 
past [13] we assume the agent correlation function C(€) to decay very fast, so that we may 
replace C{£) — > c + (5£o(l ~ c). Substitution into ([U5]) . followed by carrying out various sums 
analytically (using 8^ Ylii G{t) = x) leads us to the approximation 

2 _ n , „\ J,. „r,.\ I'' 



a 



1 

+ -(l-c) hm lim { {-^nT'-' / d/i^>(^)/i'-+'^+i Yl G{eo,h)...G{£r- 

2 t^ooV-O L^f;-^ Jo 



r,r'>0 io. 

X G(£^),/i)...G(4'_i,<0'J^/:,^^A.'^M,}L (94) 

In the second term we see that, unless r = r' = (where 5^ £i reduces directly to 
^io/'o ~ ~ ■*■)' factor Si^^i' ^ will always leave us with one residual factor 6n that is 
no longer compensated by a time summation, resulting in the removal of the corresponding 
term for 5^ 0. Hence one retains only 

1 f°° LL 1 

= 2 + ^(^^ (1 + ^^)2 + 2^'^ ~ ^) 

1 f°° u 1 

= -(l+c)| d^lQ{^l)^^^—:^ + -{l-c)Qo (95) 

For C — > 1 (strictly inconsistently fake histories), where ^>(/i) — Qo\ and Qo — > K2, 

this expression reduces to the familiar approximation formula known from inner product MGs 




Figure 5. The bid distribution P{A) = L'^ I]^<l(^[^^^(^)]) of inner product MGs 
(j78p . as measured in numerical simulations (histograms), with N = 1025 and a = 2. 
Statistics are calculated over 200. iV time steps, following equilibration. The solid lines 
(barely visible, since they are virtually coinciding with the histograms) are Gaussian 
distributions with average and variance identical to those observed, for comparison. 
Top row: f[A] = A, with C S {0, i, 1} Bottom row: f[A] = sgn(A), with C G {0, i, 1}. 
The left pictures refer to strictly real histories, the right pictures to strictly fake 
histories, with the middle ones representing a mixing of the two. These data 
(corresponding to the same experiments as those of figure [3]) are representative of the 
typical behaviour throughout the ergodic regime. The bids are always indistinguishable 
from Gaussian variables. 



(apart from various factors K2, which appear exactly as predicted by ([53]) and (j6T|) : 
5.4- Compactification of the theory and universality for C 

The final equations which describe the static limit of our theory, describing time-translation 
invariant stationary states, can be compactified further. The simple exponential shape of 
the spectrum (189p allows us to express the two relevant averages in (171|72p in terms of 
the exponential integral Ei{z) = J^dt t~^e~* [2^, resulting in the following set of coupled 
equations for the static order parameters {c, (p,x}- 

u = ' (97) 



wct[oo]^/^/J2(1 + c) (98) 
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(99) 
(100) 

(101) 



Equations (|97l98p simply define uj and u as short-hands for comphcated functions of {c, (j), x}- 
The volatihty and the factors Qci are calculated from 

1 



£7 



/i±^Ji+x(Qo-Qi)]/xQi f l+x(Qo-Qi: 



i+x(Qo-Qi) 



+ ^(l-c)Qo 
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(1 + QiX) 
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i+x(Qo-Qi) 
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c)Qo 



Qo 



/2[zy'(i-C)V2 + C252] 
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(102) 

(103) 
(104) 



All the details regarding the history definitions (viz. the sampling function f[A\ and the 
degree of 'fakeness' Q) are fully contained within the factors Qo,!- Our stationary state order 
parameter equations are now fully closed, and can be solved numerically. 

Of special interest is the limit ^ ^ 0, where the histories are strictly real. In this limit 
one can show that the persistent order parameters {(/>, c} become functions of a that no longer 
depend on the choice made for f[A\. The reason for this universality is that according to 
(|79l80p for C ^ one always has lim^^oQi = lim(;-^oQo = SDz f'^[az]. We may now define 
X = QoX and fiiid for C our order parameter equations for {c, cp, x} reducing to 

el/* Eiil/x) 



X 



X^ 



eVx Ei{l/x){2x+l)+X 



(105) 

(106) 
(107) 
(108) 

(109) 

The values of {c, (/>, x} as solved from the above closed set can only depend on the remaining 
control parameter a. Since \\m.z^Q ze^ Ei[z) = 0, the ergodicity-breaking transition point 
where x ~^ oo, corresponds as always to (p = 1 — a. We would only need to know the choice 
made for f[A] if we wished to disentangle x arid Qo from x, or to calculate the volatility which 
would be the solution of 



u = uJcj[oo]^/a/ y2{l + c) 

4> = 1 - Erf [li] 

Erf[n] =a{l-x-V/* Ei{l/x)} 
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For f[A] = A, however, we have a degenerate situation at = 0. Here we find (jllOp converting 
into 
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Figure 6. Observables c, (j) and a as functions of a = p/N , measured in numerical 
simulations for f[A\ = sgn(j4) and S = 1. Measurements are taken over 200.p time 
steps, after an equilibration of 200. p time steps; system sizes: pN"^ = (except for 
a = 16 and a = 32, where N = 1025 and N — 513, respectively). Each panel shows the 
result of four simulation runs: C = (strictly real histories, connected full circles) with 
biased and unbiased initial conditions, and C = 1 (strictly fake histories, connected 
open circles) with biased and unbiased initial conditions. Biased initial conditions: 
qi(0) = 0(1) (zero average random); unbiased initial conditions: <Zi(0) — 0(10^^) 
(zero average random). Dotted curves give the theoretical predictions for = and 
C = 1, upon assuming ergodic stationary states, with the vertical dashed line marking 
the transition value ckc below which the theory no longer applies. 



We find that (possibly apart from one degenerate case) that the only stationary solution has 
o" = 0, so also = cr^ = 0. This, in turn, tells us that x = oo for all a, so the f[A] = A 
system is always critical and the ergodicity assumption is always inconsistent. 

6. Theory versus numerical simulations of inner product MGs 

Below we test the predictions of our theory further against numerical simulations of inner 
product MGs without agents' decision noise, so T = and a[co\ = 1. Our equations allow 
for such noise, but it represents a simple generalization that is well understood and delays 
considerably the equilibration in simulations. In contrast to look-up table MGs, where the 
non-Markovian nature of the microscopic process involves delayed forces going back over a time 
interval At = 0{N-^\ogN), in inner product MGs this time interval scales as At = 0{N^). 
This already limits severely the scope of numerical simulations, and, given the multitude 
of control parameters to be varied, system sizes will have to remain modest, here mostly 
N = 0{lQ^). Even then, finite size and finite time effects remain a serious problem, especially 
for large a. 

6.1. Inner product MG with f[A\ = sgn{A) 

In an earlier section we have already compared the observed spectra of the history covariance 
matrix to our theoretical predictions, see figure HI showing quite reasonable agreement. In 
figure [6] we compare the predicted values of the observables {c, (j), a} (obtained by solving 
our order closed equations (|97m04p numerically) to numerical simulation data, for C £ {0, 1}, 
S = 1, and different values of a. We see our prediction that the transition point is independent 
of C, i.e. of whether histories are real or fake, is confirmed quite convincingly. In the ergodic 
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Figure 7. Observables c, 4> and cr as functions of measured in numerical simulations 
for f[A] — sgn(A), a = 2, and S = 1. Measurements are taken over 400. time steps, 
after an equilibration of 400. time steps; system sizes: N = 1025. Circles give the 
simulation results. Curves give the theoretical predictions, upon assuming ergodic 
stationary states. Strictly real histories and stricly fake histories correspond to C = 
and C = 1, respectively. The random deviations between theory and simulations are 
compatible with the expected finite size effects of order N^^^^ w 0.03, except for the 
volatility where fluctuations are small and the weak predicted trend is not confirmed. 



regime a > ac, where our theory apphes, we also observe a good agreement in terms of the 
values of c and </> (except for large values of a, where finite size and finite time effects become 
severe). In terms of the volatility a, for which we had to make further approximations (to 
express non-persistent terms in persistent ones) theory and experiment exhibit deviations: the 
experimentally observed value for a appears roughly independent of C in the regime where 
the data are still reliable (if not a higher value for C, = compared to C = !)> whereas 
the theory predicts a slight volatility reduction due to having real histories compared to the 
fake ones. This suggests that, at least for f[A] = sgn(^), the usual crude step to replace 
C{i) c+ (1 — c)6ui in the true formula (|93p for the volatility, resulting in the approximation 
(j95p . is less appropriate for real than for fake histories. A careful analysis of the short-time 
behaviour of the kernels C and G will be required to understand why this is so. Fortunately, 
the volatility drops out of our order parameter equations for C G {0, 1}. Upon measuring order 
parameters as a function of ("j in order to probe those cases < C < 1 where the volatility 
does play a vital role in closing our equations, one still observes excellent agreement in those 
regimes where simulation data are reliable. 

6.2. Inner product MG with f[A\ = A 

In contrast to the inner product MG with f[A] = sgn(A), for f[A] = A the theory predicts 
a strong dependence of the volatility on the parameter with a degeneracy at = 0. This 
prediction was already confirmed qualitatively by the data in figure [3l In figure [8] we inspect 
this dependence in more detail for a = 2, and find an excellent agreement between theory 
and simulations, even for the volatility (in contrast to the situation for f[A] = sgn{A), where 
the predicted volatility was somewhat lower than the observed one). According to (jllip we 
should expect that the only stationary solution for ( = will have cr = 0. In the limit 
^ (where solving our closed equations becomes nontrivial as a result of having o" — > and 
X — > oo) figure [8] indeed confirms the predicted stationary state with vanishing volatility, and 
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Figure 8. Observables c, 4> and a as functions of measured in numerical simulations 
for f[A] = A, a = 2, and 5* = 1. Measurements are taken over 400. time steps, 
after an equilibration of 400. time steps; system sizes: N — 1025. Circles give the 
simulation results. Curves give the theoretical predictions, upon assuming ergodic 
stationary states. Strictly real histories and stricly fake histories correspond to C = 
and C = 1 1 respectively. The deviations between theory and simulations are compatible 
with the expected finite size effects of order iV^i/^ ^ g.OS. 



the numerical solution of the order parameter equations confirms that x ~^ cxd. However, if 
the experiment is repeated for larger values of a one finds that beyond a critical value ~ 6 
this zero volatility stationary state is no longer stable, and the fluctuations diverge. This is 
illustrated in figure [9j In fact, this runaway solution is found to exist also for C > 0, see figure 
[Tot suggesting a C-dependent criticality a^(C)- Also at the lower end of the a-scale, below the 
conventional transition point, one finds the system generally entering the runaway state, see 
e.g. the volatility data in figure [TTl below a second critical point a~{C), provided one chooses 
unbiased initial conditions. For intermediate a values, where the volatility remains finite, our 
theory is found to predict the volatility a correctly (within the accuracy limitations imposed 
by finite size effects) . In the right panel of figure [11] we show the result of estimating from 
numerical simulations the location of a^(C) in the (a,C) plane, which appears to approach 
= 1 for a — > 0; determining this line experimentally is not entirely trivial in view of the 
(expected) nontrivial bid fluctuations close to this line, see eg. figure [12j Carrying out a 
similar exercise for the location of the instability a^(C) in the high a region is unfortunately 
ruled out, due to the extreme relation times required for large a. In summary: for sufficiently 
small C (i.e. sufficiently 'real' histories), the usual high volatility state of conventional MGs 
in the non-ergodic regime becomes a runaway solution in inner product MGs with linearly 
sampled histories (a true instability with infinite volatility), which for C = vanishes slightly 
below the standard MG's x = oo transition a = Oc if one increases a, but then re-appears at 
a larger value a^iC) and there becomes the only solution. It must in principle be possible to 
calculate the critical values a^(C) from our order parameter equations; in practice, however, 
this requires solving our equations for the correlation and response functions (formulated in 
terms of both the effective agent process and the effective overall bid process) at finite times, 
which at present we are unable to do. 
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Figure 9. Overall bids A{(.) versus time t — 1 . . . 300,000, as measured in numerical 
simulations of the inner product MG with f[A\ = A and C = (fully real histories), in 
systems of size N — 1025. Top row: unbiased initial conditions, qi{0) — 0(10"^) (zero 
average random); bottom row: biased initial conditions, qi(0) = 0(1) (zero average 
random) . The data indicate that there exists a critical value a J such that for a > aji" , 
rather then evolve towards tr = 0, the system is unstable and the bid fluctuations 
diverge. If the (weak) dependence on initial conditions in the graphs for a = 6 is not 
a finite size effect, then there is a parameter regime where the zero volatility state and 
the 'runaway' solution coexist. 

7. Discussion 

In this paper we have generalized the theory of look-up table Minority Games with real market 
history [11] to a larger family of models, and applied it to MGs with so-called inner product 
strategy definitions. The latter MG versions had, surprisingly, never been solved; not even 
in their simplest Markovian version where the histories are fake. At a mathematical level it 
was not a priori clear which form the inner-product MG theory would take: in [11] it was 
found that the key object in the theory of look-up table MGs with real market history was the 
so-called history frequency distribution, but already on simple scaling grounds it is clear that 
this object cannot be defined for inner product MGs. It is satisfactory to find in the present 
study how the generating functional formalism resolves the issue: a more general quantity 
takes over the role of the history frequency distribution, viz. the eigenvalue spectrum of the 
history covariance matrix, which reduces to the former for look-up table MGs but is also well- 
defined for inner product MGs. This resolution involves a generalization of the short history 
correlation time ansatz that led in [11] to closed stationary state order parameter equations, 
resulting (after a random matrix spectrum calculation) also for the present inner-product MGs 
in closed equations for observables such as the persistent correlations, the fraction of 'frozen' 
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Figure 10. Overall bids A{£) versus time £ = 1 . . . 300,000, as measured in numerical 
simulations of the inner product MG with f[A] = A, in systems of size N = 1025 
and following unbiased initial conditions. Here C — 0.1. The unstable state with 
diverging volatility exists also for > 0, with the instability setting in at a critical 
value a < that increases with (. Also the characteristic time scale for the 

instability to manifest itself appears to increase with C (note the different time scales 
of figures [Hand [To]). 

agents, the integrated response, and the market volatility (although the formula for the latter 
as always involves further approximations). 

We find, provided all the relevant integrals and averages exists, that the phase diagrams 
of the two MG versions are identical. However, we encounter interesting differences in terms of 
the models' static and dynamic phenomenology, dependent upon which functions f[A] of the 
past overall market bids A are being sampled in the inner product MGs. For bounded functions 
such as f[A] = sgn(^) one continues to find behaviour that, although quantitatively different, 
is qualitatively similar to that of look-up table MGs. If the histories are strictly real, the 
theory even predicts complete universality of the static observables, independent of the choice 
made for f[A]. In contrast, for unbounded sampling functions such as f[A] = A the situation 
is quite different. There can now be a profound dependence of the volatility on the degree to 
which the histories are real. Especially when the histories are strictly real and f[A] = A we 
find that the inner product MG has only one stationary solution, where the volatility a is zero 
and the model is structurally critical (i.e. x = oo). In terms of the usual control parameter 
a = p/N, we find that in the regime a > Oc (the standard non-ergodicity transition of the 
MG) this state (cr, x) = (0, cxd) can only be reached for a < a'^{C), with a^(0) 6 and a^(C) 
increasing with For a > a'^iC) the system destabilizes and the fluctuations diverge. In the 
low Q regime one finds a similar critical line a~ (C) marking a destabilization transition as one 
reduces a further, with a = oo for q < a~{C), but here this destabilization occurs only for 
unbiased initial conditions. 

Those technical questions that remain open in the present study are the familiar ones that 
are similarly answered for all previous MG models. They relate mostly to our inability so far 
to extract from our order parameter equations exact solutions for time-dependent observables 
(beyond just a few time steps), or exact stationary state solutions in the non-ergodic regime. 
In addition one would like to explore further the possibility of finding solutions for MG versions 
in which the effects of having real histories are more dramatic; both the present study and 
its predecessor rely on expansions that require history correlations to be weak. Finally, 
we would like to point out that the observed instability of the f[A] = A inner product MG 
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Figure 11. Left: volatility as a function of a = p/N, measured in numerical 
simulations for f[A] = A and 5=1. Measurements are taken over 200. p time steps, 
after 200.p equilibration time steps; system sizes: pN"^ = 2'^" (except for a = 16 and 
a = 32, where A*" = 1025 and N = 513, respectively). We show the results of four 
runs: C = 0-1 (predominantly real histories, connected full circles) with biased and 
unbiased initial conditions, and C = 1 (strictly fake histories, connected open circles) 
with biased and unbiased initial conditions; initial conditions are as in previous figures. 
Dotted curves give the theoretical predictions for ^ = 0.1 and C = 1; upon assuming 
ergodicity. The system destabilizes {a = oo) for a > Q!^(C) (any initial conditions, 
terminating branches at a « 9) and for a < a~ (unbiased initial conditions only, 
terminating branch at a ~ 0.2). Right: estimated critical line a~{() in the low a 
region of the (a, () plane, marking the onset of the ct = oo instability, obtained by 
numerical simulations with N = 2049 (C < 0.77) and N = 4097 (C > 0.77). The 
marker sizes indicate the error bars. The vertical dashed lines in both figures mark 
the X ^ 00 transition value ac below which the ergodic theory no longer applies. 

is qualitatively different from what has been observed in previous MG versions. Although in 
its present form this specific model is ill-constructed, it could perhaps form the basis of new 
models aimed at explaining more robustly the so-called stylized facts in real markets, since 
its structural instability is not a finite size effect and therefore requires no careful tuning of 
control parameters. The phenomenology observed for f[A] = A is novel, should be expected 
to emerge also for other unbounded choices for the sampling function f[A], and emphasizes 
even more the need for progress in solving the GFA order parameter equations for finite times. 
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